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MODELING OF TRANSIENT HEAT PIPE OPERATION 

1. Introduction 

The use of heat pipes is being considered as a means of reducing the peak 
temperature and large thermal gradients at the leading edges of reentry vehicles 
and hypersonic aircraft, and in nuclear reactors [1,2]. In the basic cooling concept, 
the heat pipe, which is a highly effective heat transport device, covers the leading 
edge, a portion of the lower wing surface, and a portion of the upper wing surface. 
Aerodynamic heat is mainly absorbed at the leading edge, and transported through 
the heat pipe to the upper and lower wing surface, where it is rejected by thermal 
radiation and convection. Once fully operational, the near isothermal heat pipe 
virtually eliminates temperature gradients and reduces peak temperatures. 

A previous feasibility study[l] of heat pipes for this application recommends a 
rectangular cross-section for the heat pipe based on weight per unit surface area of 
heat pipe cooling structure and fabrication considerations. A schematic diagram of 
the physical model based on results presented in referenece[l] is shown in Figure 1. 

Previous experimental observations [2-8] suggests the following sequence of ev- 
ents during heat pipe startup from the frozen state. Initially the working fluid is 
in the solid state and the vapor density is extremely low so that free molecular 
flow conditions prevail throughout the vapor space. The input heat flux over the 
evaporator starts to melt the frozen substance in this region, and the heat transport 
from the hot zone to the adjacent zone proceeds quite slowly via axial conduction 
through working fluid and capillary structure, and heat transfer in the vapor is 
almost negligible. 

When energy is continuously added to the evaporator, the frozen working fluid 
in the evaporator is melted so that evaporation can take place at the liquid-vapor 
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interface and the vapor can flow into the condenser section due to the large pressure 
gradient. Hence, vapor freezes on the inner surface of the frozen working fluid in the 
cold zone so that the vapor-solid interface temperature increases until the melting 
temperature is reached. During this stage, energy is mainly transferred as latent 
heat owing to vaporization in the heated zone, and condensation and freezing in the 
cooled zone. Also, vapor flow may be choked at the exit of the evaporator because 
of very low pressure in the cold zone. 

This process continues until the frozen working fluid is completely melted and 
the continuum flow regime reaches the end of the heat pipe, at which time liquid 
returned to the evaporator is sufficient for normal transient operation. Eventually 
the heat pipe may reach a steady state condition. As suggested, during the startup 
of the heat pipe from a frozen state, the behavior of vapor flow may be divided into 
three distinct phases for convience of analysis. 

Phase I: The vapor flow in the heat pipe is in free molecular flow through the 

vapor space. 

Phase II: In the vapor space, a region of continuum flow is established in the 

heated zone and a continuum flow front moves toward the heat pipe 
end. The vapor flow may be choked at the end of the evaporator. 
Phase III: Continuum flow exists over the entire heat pipe length in the vapor 
region and the sonic limit is not encountered. 

2. Mathematical Model Development 

On the basis of experimental observations, basic governing equations axe writ- 
ten to determine the startup, transient, and steady state performance of a heat pipe 
which has initially frozen alkali-metal as the working fluid. These equations can be 
coupled by several types of boundary conditions on the heat pipe surface, such 
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as specified temperature, heat flux, convection and radiation boundary conditions. 
The boundary condition at the liquid-vapor interface depends on the three phases 
of vapor flow dynamics mentioned in section 1. 

2.1 Transition temperature 

Continuum flow in the vapor space is considered to be established when the 
mean free path, A, is substantially less than the minimum dimension,!), of the vapor 
flow passage, e.g., 

K,= ^< 0.01 ( 1 ) 

From kinetic theory of gases [9], the dynamic viscosity, and the mean molecular 
velocity can be expressed as 

n = 0A99p\V (2) 


v = 



( 3 ) 


Elimin ating the mean free path from equations(2) and (3) yields the tempera- 
ture of the vapor space corresponding to the given mean free path, 


T* > 


i r M ( p \ 2 

2x 10 ~ a Y u \Jd) 


( 4 ) 


where p is the density of the vapor, p is the dynamic viscosity of the vapor, R u is 
the universal gas constant, and M is the molecule weight. 

Iterations are required to obtain a value of T* due to the dependence of temper- 
ature of properties. Figure 2 illustrates the transition temperature, T*, from free 
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molecular to continuum flow as a function of minimum vapor passage for sodium. 
When the temperature of the vapor space is greater than that calculated by Equa- 
tion^), continuum flow is assumed to be established in the vapor space. 

2.2 Heat pipe shell and capillary structure with fluid 

Unlike those for a conventional cylindrical heat pipe, aerodynamic heating will 
cause the highest temperature to occur at the outer shell of heat pipe. Thus, to 
simplify the analysis, it is assumed that the inner shell is sufficiently thin to ne- 
glect its thermal resistance and capacitance. Furthermore, the following additional 
assumptions are made. The thickness of the heat pipe and the rib is assumed to 
be much smaller, than the width, and temperature gradients would be developed 
primarily in the chordwise direction during heat pipe startup. Therefore, the rib 
may be neglected in the analysis. Also, the heat pipe is assumed to be symmetric 
about the stagnation line so that the upper half section of the heat pipe is consid- 
ered. Hence, a two-dimensional model is considered for mathematical formulation, 
and a simplified cross section is shown in Figure 3. 

The unsteady, two-dimensional conduction equation is applied to the heat pipe 
shell and the capillary structure to evaluate the temperature distributions under 
the following assumptions: 

1. The heat transfered through the wick and working fluid is by conduction 
only; 

2. The presence of the wick structure does not affect the melting of the work- 
ing substance, and melting occurs over a small, finite temperature range, 
AT, around the melting temperature; 

3. The boundary conditions at the liquid-vapor interface are subject to the 
phases of vapor flow dynamics noted in section 1. 




Figure 3. Simplified cross-section A-A of the heat pipe cooled leading edge. 
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4. The radiation heat transfer at the liquid-vapor interface is neglected based 
on the small thermal emissivity of liquid sodium measured [10]. 

The governing equations and boundary conditions for the heat pipe shell and 
capillary structure with working fluid are expressed in the same form: 



Ti = T 0 for t < 0 (6) 

~ K i-q£- — Q(X,t) at Y = f 1 (X) for 0 < < 4> e (7) 


Kl di ~ kc 


•(Ti - T cr ) + <t£(Ti* - T r 4 ) at 


Y = fi(X) for ij> e < 


( 8 ) 


T T A V ^ If 

T\ = Ti and Ai— = K 2 -z— 
On on 


at Y = fiiX) for 0 < ^ < ^> c (9) 


dTi 

dS 


= 0 


at rp = 0 and rj; = ij> e 


( 10 ) 


where the subscripts 1 and 2 denote the heat pipe shell and the capillary structure, 
n represents the unit outward normal direction, and S represents the chordwise 
direction. 
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In the X - Y coordinates 

( 11 ) 

where l x and l y are the direction cosines between the surface outward normal, n, 
and the X, and Y axes, respectively. 

2.3 Analysis of the vapor flow 

Analysis of the vapor flow is necessary to provide the boundary condition at 
the liquid-vapor interface when continuum flow is established in part of the vapor 
space. However, the behavior of the vapor flow as mentioned in section 1 is very 
complicated due to the extremely small vapor density of the metallic working fluid 
at low vapor pressure and the large pressure gradient in the axial direction. Limits 
on the vapor flow are encountered and considerable thermal resistance exists at 
the liquid-vapor interface due to phase change of the working fluid. The case of a 
flat plate heat pipe with asymmetrical boundary conditions shown in Figure 3 is 
considered for the vapor flow passage. 

Evaporation and condensation 

From kinetic theory, a significant thermal resistance exists at the liquid-vapor 
interface for liquid metal[ll] and this resistance increases with decreasing vapor 
pressure. This means that the difference between the vapor temperature and the 
interface temperature is not negligible. At the interface, the interphase mass transfer 
has been stated from the viewpoint of kinetic theory as a net molecular flux which is 
the difference between the rate of arrival of molecules from the vapor space towards 
the interface and vice-versa. While condensation is proceeding, the arrival rate of 
molecules exceeds the departure rate. During evaporation, the reverse takes place, 


.dT 

dn 


dT 

ax' 


dT 


K — — K -f K-^yty 
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and during the equilibrium state, a net molecular flux is zero. Hence, evaporation 
and condensation are modeled by using an expression for the net rate derived from 
the kinetic theory of gases[12]: 




/ 2 a \ I M 

~ \2 — a/ V 2 irR u 


ly/T f y/T v \ 


(12) 


where a is the condensation or evaporation coefficient which is assumed to be unity, 
rh 0 is the rate of condensation or evaporation per unit area, M is the molecular 
weight, R u is the universal gas constant, P/ and T/ are pressure and temperature, 
respectively, at the interface, and P v and T v axe pressure and temperature of vapor, 
respectively. 

Limitations of vapor flow 

After continuum flow is established in the vapor space, because of the low 
density of vapor at low pressure and the large pressure gradient, the vapor flow 
is choked at the end of the evaporator, even for a low heat flux, just as it is at 
the throat of a convergent nozzle. This sonic limit is the first amoung several 
limitations encountered, and the performance of the heat pipe is restricted until the 
vapor temperature increases significantly so that the velocity of the vapor leaving 
the evaporator is less than the sonic velocity. 

The expression for the sonic limit[6,13] in terms of flow conditions at the begin- 
ning of the evaporator by energy balance on a control volume enclosing the entire 
evaporator is: 

= po Ah fgy fiR u T o 
v/2( 7 + 1 )M 


where A is the cross-section area of vapor space, hf g is the enthalpy of vaporization, 
7 is the ratio of specific heat, M is the weight of a molecule, iZ u is the universal gas 
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constant, T a is the vapor temperature the beginning of the evaporator and 70 is 
taken as the saturation density corresponding to T 0 . Experimental verification of 
Equation(13) has been made by Kemme[14] for sodium, potassium and cesium heat 
pipes. 

Figure 4 shows the sonic limit, the entrainment limit [15] and the Reynolds 
Number for a heat pipe which has a cross-section area of 0.55 cm 2 , hydraulic radius 
of .32 cm, 100 mesh screen wick, and sodium as the working fluid. The type of 
limitation restricting performance of heat pipe is determined by which limitation 
has the lowest value at the temperature under consideration. Thus, the first limit 
encountered is the sonic limit, as shown in Figure 4. However, when the actual 
axial heat transfer required is below these limits, no limits axe encountered. 

Modeling of vapor flow during phase II 

Even though some of the working substance is in the liquid state, the transition 
temperature from free molecular flow to continuum flow is much greater than the 
melting temperature, so the vapor is assumed to be free molecular during phase I, 
and an adiabatic boundary condition is applicable at the interface. However, during 
phase II, a region of continuum flow is assumed to be established in. the adjacent 
vapor space when the interface temperature is greater than the transition tempera- 
ture, while in the cold zone the vapor is still in free molecular flow. The continuum 
flow region grows until it reaches the end of the heat pipe. The temperature in the 
region of free molecular flow remains unchanged except in the vicinity of the con- 
tinuum flow region. Therefore, an imaginary plane, which is adiabatic and normal 
to the liquid-vapor interface, is assumed to divide the two vapor flow regions at the 
point of the transition temperature, and the dividing plane moves toward the end of 
the heat pipe as the location of the transition temperature at the interface moves. 


Axial Heat Transfer Limit x 10, w/cm 

10.0 15.0 20.0 25.0 30.0 35.0 40.0 45.0 50.0 
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Figure 4. Heat pipe vapor flow limits and correspondinf Reynolds Number. 


Reynolds Number x 10 
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In the continuum flow region, energy is mainly transported as latent heat of 
the working fluid, and variations of temperature and pressure in the axial direction 
are significant. Even though continuum flow exists in the vapor space, the vapor 
pressure is low and the pressure gradient in the axial direction is large, so the heat 
transfer is limited by the choked flow condition, and supersonic vapor flow and a 
shock front[6,7,16] may occur in the condenser. 

In this research, the overall performance of a heat pipe is of more interest than 
the details of one portion, and as noted, the ultimate heat transfer rate in the axial 
direction is limited by the sonic limit. Thus, the vapor flow during phase II may 
be approximated without loosing the generality and accuracy by evaluating the 
sonic limit properly. In order to evaluate the sonic limit, the total heat input to 
the vapor space, which can be obtained by applying Equation(12) to the element 
adjacent to the continuum flow region, is equated to the sonic limited transport 
from Equation(13) as follows: 


^ / M Pf i P v w _ P A c y 7 R u T v 

*K{y/Ti' i y/Kl - V2(t + 1 )M 


(14) 


where A c is the cross-section area and W is- the width of the vapor space. From 
Equation(14), the uniform vapor temperature can be computed by iteration. This 
vapor temperature may be lower than true vapor temperature at the beginning of 
the evaporator, but at this stage the variation of density with temperature and the 
density itself is very small, so the variation of sonic limit with respect to temperature 
is also small, as shown in Figure 4. Thus, the heat flux at the liquid-vapor interface 
and the sonic limit can be obtained by solving Equations(12) and (13), respectively, 
with the vapor temperature obtained from Equation(14). This method may be used 
until the vapor flow state reaches phase III. 
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Vapor flow during phase III 

In this stage, the entire working fluid is completely melted and continuum flow 
exists throughout the vapor space. However, the heat pipe has not reached the 
desired operating condition and the density of the vapor is still small, so that com- 
pressibility should be considered. Also the amount of energy stored in the vapor 
space is negligible. Thus, quasi-steady, compressible, one-dimensional laminar flow 
in the vapor space is considered. For purposes of formulating the mass, momen- 
tum, and energy equations in one-dimensional form, the velocity is taken to be 
the average velocity, which is approximated by the velocity distributions based on 
the similarly solutions of semiporous channels. In addition, friction at the liquid- 
vapor interface, variations of vapor quality, and momentum and energy factors are 
similarly calculated. The vapor is evaporated at the interface with mass injection 
rate,m 0 , per unit area. It is assumed that this vapor flows inward with a normal 
component of velocity only, and joins the axial flowing stream. 

The principles of conservation of mass, momentum, and energy in the axial 
direction, and the assumptions noted, yield: 

D 7^(P V ) = Th o (15) 


dP d . o \ FpV 2 

dZ + dZ ^ ap ^ ~ 8 D 


( 16 ) 


D± 

dZ 


P V(K + ^) 


= m 0 (h o + V 0 2 ) 


( 17 ) 


where D is the height of the vapor space, h v is the vapor enthalpy and h 0 is the 
vapor enthalpy at the interface. 
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Friction factor F of the surface is written as 


F = 


8 r v 

pV 2 


(18) 


Momentum and energy factor, a and /?, respectively, are expressed as 

D 

“ = ok/ ^ (19) 

0 

D 

P = v ' d v < 2 °) 

0 

Normal velocity of the vapor at the interface is expressed in terms of heat flux 
and latent heat of vaporization as follows: 


F 0 = 


Q 

hfgPoA 0 


( 21 ) 


where Q is the heat input rate at the interface, p Q is the density for the interface 
temperature and A a is the interface area. 


Coupling vapor flow effects 

When the Mach number is less than about 0.2, compressibility is neglected 
and friction effects at the interface may be negligible due to the low velocity of the 
vapor. Thus, most studies[17,18] of heat pipe performance have assumed that the 
temperature is uniform throughout the vapor space. This approximation eliminates 
difficulties encountered from solving the vapor flow dynamics and gives simple re- 
sults for the case of a low operating temperature with the small heat flux. When 
the vapor temperature is uniform in the vapor space, no thermal resistance exists, 
so a maximum amount of energy can be transfered through the vapor space for the 
given operating condition. 
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The energy stored in the vapor space is negligible due to low density, so the 
energy entering the evaporator is equal to the that leaving the condenser. Then 
Equation(12), which describes evaporation and condensation, is applied to every 
element at the interface, and the energy balance in the vapor space is: 


£ 




p, ' 

vz. 


( 22 ) 


where m e is the number of elements at the interface in the evaporator and m c is the 
number of elements in the condenser. Equation (22) can be solved by using iterations 
to obtain uniform vapor temperature as long as the temperature at the interface is 
known. Also, the heat flux at the interface is evaluated using Equation(12) with 
known temperatures. 

However, when liquid metal is used as the working fluid and the heat flux ap- 
plied to the surface is large, the velocity of the vapor is large, so the temperature 
drop in the vapor space should be considered. For a heat pipe which has a rect- 
angular cross-section area of 1.77 cm 2 , Figure 5 shows considerable temperature 
drops for large heat fluxes. The temperature drops in the vapor region imply that 
the thermal resistance causes less energy to be transferred from the evaporator to 
the condenser. Thus, the evaporator temperature is higher and the temperature of 
condenser is lower than in the other case. 

When the vapor flow dynamics are coupled with the heat conduction in the heat 
pipe shell and capillary structure at the interface, the governing equations are solved 
simultaneously -with unknown boundary conditions at the interface, so iterations are 
required for every time step until both results match at the interface. This method 
may yield accurate results, but in general consumes much computational time due 
to the iterations. Also convergence may not be reached. Hence, an approximation 
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method is employed to eliminate these difficulties. 

Instead of simultaneously solving governing equations for both regions, govern- 
ing equations for the vapor region are separately solved with given a heat flux, so 
that the temperature drop, AT,,, can be obtained. The thermal resistance,#,,, in 
the vapor space may then be evaluated from Equation(22): 

AT 

Q=~iT ( 23 ) 

In order to get thermal resistance, the total heat flow rate is needed, but without 
solving both governing equations, the heat flow rate is also unknown. Thus, the 
known heat flow rate at previous time step is used, and heat flow rate evaluated at 
the present time step would be used at the next time step. 

When this evaluated thermal resistance is imposed at a different place without 
changing total heat transfer rate, the uniform vapor temperature can be used in the 
vapor space. Since the total heat transfer rate through the vapor space, which is 
computed by using Equation(12), is the same as the heat input at the liquid-vapor 
interface in the evaporator or heat output at the interface in condenser, a certain 
fictitious layer which has the same thermal resistance as that evaluated in the vapor 
space may be placed at the liquid-vapor interface. Thus, the temperature at the 
new interface is lower than that at the true interface due to the resistance in the 
fictitious layer for the evaporator, and vice-versa for the condenser. 

In order to compute the uniform vapor temperature, the temperatures at the 
new interface axe used for Tf i in Equation(22), and the heat flux at the interface 
is calculated by using Equation(12) with the uniform vapor temperature and the 
temperatures at the new interfaces. This heat flux then serves as the boundary 
condition at the liquid-vapor interface when solving the governing equations for the 
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heat pipe shell, capillary structure, and vapor flow at the next time step. By using 
this approximation scheme, iterations can be excluded while the vapor flow effects 
are retained. 


3. Results and Discussion 

To verify the model developed for heat pipe cooled leading edges, numerical 
results axe compared to Camarda’s experimental results[2]. Since the shape of the 
cross-section of Camarda’s model is different from that of the present model, the 
rectangular cross-section for this research is modeled such that the cross-section ar- 
eas occupied by the heat pipe shell, the capillary structure, and the vapor space are 
the same as for his experiment. Also, the surface area where the heat flux is applied 
is chosen identical to that of the experiment so that the width and length of heat 
pipe are identical, but the thickness of the component between the two heat pipes 
is different. However, the thickness is much smaller than other dimensions, and the 
temperature drop at the cross-section is negligibly small compared to the temper- 
ature drop in the axial direction, so the effects of the difference in the thicknesses 
may be minimal. Thus, effects of energy storage and fusion of working fluid in the 
heat pipe shell and the capillary structure are the same, and the effect of vapor flow 
dynamics may be simulated. Since the minimum dimension of the vapor space is 
different due to a different thickness of the heat pipe, the transition temperature 
which is evaluated based on the dimension of the vapor space is not identical. Also, 
this temperature determines the boundary condition at the liquid-vapor interface, 
so to simulate the Camarda case, the inner diameter of the circular heat pipe is 
used to calculate the transition temperature. 

The material of each component used in the computation is the same as in the 
experiment except that the Braze alloy is assumed to be hastelloy x to simplify the 
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heat pipe shell. The properties of the materials are varied with temperature during 
numerical calculation. 

Heating distribution on the heat pipe surface is modeled based on the normal- 
ized heating distribution [2] shown in Figure 6, and time stepped heat input at the 
stagnation point. The stepped heating distribution due to aerodynamic heating is 
shown in Figue 7. In addition to aerodynamic heating, a radiation boundary con- 
dition is used to extract energy from the entire heat pipe surface so that the net 
rate affects heat pipe internal operation. The section which has a positive net rate 
is considered the evaporator, and the section with a the negative net rate is the 
condenser. 

Figure 8 shows the two-dimensional grid system, dimensions, boundary con- 
ditions, and materials used to represent a leading edge cooling heat pipe. It is 
assumed that the temperature of the heat pipe is initially the ambient temperature, 
which is below the melting point of the sodium working fluid. Also, a transition 
temperature of 700 °K is used. 

An explicit method with a few implicit steps is used as the time stepping 
scheme. Until continuum flow is established in the vapor space, a constant time 
step of 10 seconds is used. This time step is gradually decreased due to increasing 
the stepped heat input on the external surface. 

Figure 9 shows a comparision of the numerical and experimental results along 
the heat pipe for different times. Detailed numerical results are plotted on Figures 
10 and 11. Since most of the energy, which is transferred by conduction through 
heat pipe shell and capillary structure in the axial direction, is added near the 
stagnation point, the temperature near the stagnation point increases rapidly, and 
a large temperature gradient is established in the heat pipe. This implies that the 
heat pipe is not effective due to the extremely small vapor density up to about 550 
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seconds. Temperatures in the rest of the heat pipe are not much different from the 
initial condition. At 550 seconds, the temperature at the stagnation point is greater 
than the transition temperature, and temperature gradient at the region adjacent 
to this point decreases. This means that a continuum flow region is established 
in the vapor space, so energy is mainly transported by the latent heat. However, 
most of the working fluid is in the frozen state and free molecule flow prevails in 
the vapor space. As time increases, the continuum flow region is expanded in the 
vapor space so that at about 750 seconds one-third of the heat pipe performs as a 
heat pipe. The rest of the heat pipe is still like a solid in which heat is transported 
only by conduction. Between the two regions, a large temperature gradient exists, 
and the continuum flow front is located here. At 1000 seconds, most of the heat 
pipe is active except near the end of the condenser section, and the working fluid 
is completely melted. After 1500 seconds, the heat pipe reaches a fully operational 
condition. 

Figures 10 and 11 show temperature distributions at the external surface and 
the interface, respectively, obtained by numerical analysis. Figure 12 shows the 
vapor temperature variation during startup. The vapor temperature increases cor- 
responding to increasing heat flux. However, the vapor temperature increases very 
slowly because most of the heat added in the vapor space of the evaporator is 
extracted at the interface of the condenser to heat up the adjacent cold zone. 

Even though the numerical results do not match exactly the experimental re- 
sults due to different configuration, material properties and boundary conditions, 
it is clear that the model does predict approximately the correct startup-time and 
approximately the correct temperature distribution. It is important to note that 
heat inputs to the model were based on aerodynamic heating computation which 
were only approximated in the experiments. 


28 

Finally, to show the effectiveness of the heat pipe for cooling a leading edge, 
numerical calculations were executed for the same example problem, except that 
an adiabatic boundary condition was applied to the interface during the entire 
startup. Development of large axial temperature gradients is identical due to the 
free molecular condition in the entire vapor space in the first part of transient heat- 
ing. However, beyond this stage, peak temperatures are reduced and temperature 
gradients disappear when the heat pipe is active. Figure 13 shows temperature 
distributions for both cases at a time of 1100 seconds. It is clear that the heat pipe 
is very effective in reducing temperature near the leading edge. 




Figure 12. Vapor temperature variation during startup 
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LIST OF SYMBOLS 


A 

A c 

C 

D 

F 

h 

her 

h fg 

/<(*) 

K 

K n 

L 

l 

M 

m 

m 

n 

P 

Q 

Q 

R 

R u 

S 


condensation or evaporation coefficient 
area 

cross-section area of the vapor space 
volumetric heat capacity 
height of the vapor space 
friction factor 
enthalpy 

heat transfer coefficient 
latent heat of vaporization 

functions which describe the curvatures of the wing section 

thermal conductivity 

Knudsen number 

total length of heat pipe 

direction cosine 

molecular weight 

rate of condensation or evaporation 

number of element at the interface 

unit outward direction 

pressure 

heat input rate 

heat input rate per unit area 

thermal resistance 

universial gas constant 

chordwise direction 


34 


t 

time 

T 

temperature 

T 0 

initial temperature 

U 

x - component of velocity 

V 

y - component of velocity 

X 

coordinate direction 

Y 

coordinate direction 

w 

width of the vapor space 

Greek symbols 


a 

momentum factor 

0 

energy factor 

7 

ratio of specific heat 

£ 

emissivity 

A 

length of mean free path 

P 

dynamic viscosity 

P 

density 

<7 

Stefan-Boltzmann constant 


angle used to identify chordwise 

T 

shear stress 

Subscripts 


1 

heat pipe shell 

2 

capillary structure 

c 

condenser 

cr 

convection 

e 

evaporator 

f 

liquid state 


1 


index 


o 

r 

x 

y 

V 

Superscripts 

* 


liquid-vapor interface 
radiation 
x direction 
y direction 
vapor state 


transition from free molecular flow to continuum 


